Granular Matter manuscript No. 

(will be inserted by the editor) 



Motion of grains in a vibrated U-tube 



J.R. Darias • I. Sanchez • G. Gutierrez 



o 



Received: date / Accepted: date 



o 



o 
o 



> 
in 

<N 

o 
o 



Abstract We investigate experimentally the behavior 
of the rate of growth of a column of grains, in a partially 
filled vertically shaken U-tube. For the set of frequencies 
used we identify three qualitatively different behaviors 
for the growth rate 7 as a function of the dimensionless 
acceleration F: 1) an interval of zero growth for low F 
with a smooth change to nonzero growth, analogous to a 
continuous phase transition; 2) a sigmoidal region for 7 
at intermediate values of the dimensionless acceleration 
F; and 3) an abrupt change from high values of 7 to zero 
growth at high values of F^ similar to a first order phase 
transition. We obtain that our data is well described 
by a simple differential equation for the change of the 
growth rate with the dimensionless acceleration of the 
vertical vibrations. 

Keywords Granular material • U-tube • Vertical 
vibration • Transport • Instability. 
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like heaping and granular transport in vibrated gran- 
ular deep beds in partitioned containers has been the 
subject of research for decades [1,2 ,3 ,4 ,5 ,7 ,8 ,^. In this 
work we concentrate in the experimental characteriza- 
tion of the rate of growth of the rising column of grains 
in a partially filled U-tube for different frequencies. The 
mechanism responsible for this interesting behavior is 
still an open issue. 

In a previous work [5 the behavior of a vertically 
vibrated granular system in a partially filled U-tube at 
low frequency was studied. A collective granular trans- 
port upward through one of the branches of the tube 
was observed for small grains. The experimental results 
were compared with a model based on the idea of cyclic 
fiuidization [6 . In that model it was assumed that for 
low frequencies and sufficiently high amplitudes of os- 
cillations the granular bed ffuidizes cyclically in such 
a way that an effective upward acceleration acts on 
the grains while these are in a ffuidized state. Con- 
sequently, an instability appears and the free surface 
of the bed on either branch of the U-tube rises while 
the other goes down. Therefore, one of the branches of 
the tube is rapidly filled and simultaneously the other 
empties. That model captures some relevant aspects of 
the observed behavior for low frequencies, as for exam- 
ple the exponential growth of the granular column and 
the monotonic increase of the growth rate with the di- 
mensionless acceleration of the vibration; however, we 
have found some important quantitative discrepancies. 
The above model for the measured growth rate 7 ver- 
sus the maximum dimensionless acceleration F of the 
container provides only a crude approximation for low 
frequencies. As the frequency is increased above 10 Hz^ 
it becomes clear that the model is insufficient to de- 
scribe this phenomenon. 
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In this paper we extend the work in reference [5 by 
exploring the growth rate of the height difference be- 
tween the two free surfaces, for a wider range of frequen- 
cies. We obtain a differential equation for the change of 
the measured growth rates with the dimensionless ac- 
celeration of the vibrations. 



2 Experimental setup 

For the vibrational system we used a function genera- 
tor (GW Instek SFG-2110) linked to an audio ampli- 
fier (Crown XLS 202) feeding two coupled 1000 Watts 
loudspeakers. For the granular material we used glass 
spheres with diameter between 250 and 300 fim and 
bulk static density pg = 1440 kg /m^. The top of the 
container was made permeable to air but not to grains 
(using a 125 jam nylon mesh). The container is a U- 
shaped tube 200 mm tall, with square internal cross 
section of 400 mm^. The total mass of grains used was 
124 ^, chosen to completely fill one of the U-tube branches, 




Fig. 1 Snapshot sequence of the partially filled U-tube in a typ- 
ical experiment. The grains migrate to the branch of the tube 
with the initially taller granular column. The images shown cor- 
respond to an experiment for which the U-tube was shaken at 
a frequency of 20 Hz and F = 4.5. The increase of the growth 
velocity with time can be readily appreciated. 



Images of the experiments were recorded using a 
digital camera (Pixelink PL-B741F). Filming against a 
back light (see Fig.[T]), allowed us to analyze the images 
using a commercial software to obtain the difference in 
height Ah between the two free surfaces as a function 
of time. 

The excitation signal used was a sine wave. The 
maximum amplitude of oscillation that could be reached 
was (11.1 ±0.1) mm. The excitation was controlled con- 
sidering the maximum dimensionless acceleration de- 
fined diS r = A(27r/)^/^, where g denotes the accel- 
eration of gravity, A the amplitude of the oscillations 
and / the frequency of the oscillations. The value of F 
was varied from 1 to 5.75. Lower values did not gen- 
erate any instability and higher values of this quantity 
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Fig. 2 Typical time dependence of the height difference Ah, 
between the free surfaces at the two branches of the tube. This 
plot is for f = 15Hz and T = 2.75. 



were not possible to achieve due to the limitations of 
the equipment. 



3 Results 

The height difference Ah between the granular level in 
both branches of the U-tube grows exponentially with 
time (see Fig.[2j). The growth rate 7 is determined from 
the fit of the exponential Ah = Ahoe'^^^^ where Aho is 
the initial height difference and t is the time. 

The graph of the growth rate 7 versus dimension- 
less r is shown in FiglS] The horizontal axis begins at 
r = 1^ because for < 1, 7 = 0. For each curve the os- 
cillation frequency is fixed and the amplitude is varied. 
For all frequencies investigated there was a seemingly 
continuous transition from zero growth at the lowest 
amplitudes to a finite 7 at higher amplitudes. It is ap- 
parent that for the larger frequencies (15 Hz and 20 Hz) 
the growth rate reaches a saturation value, for suffi- 
ciently large amplitudes. For three of the intermediate 
frequencies {10 Hz^ 12.5 Hz and 15 Hz) we see that 7 
reaches a point where it changes abruptly to zero at a 
threshold of F (see arrows in Fig. [3]). This suggests the 
occurrence of a first order phase transition, but more 
work needs to be done to characterize this observed 
abrupt transition. For the lowest frequency measured 
(7.5 Hz) we could not determine whether there is a sat- 
uration region or a transition from high nonzero growth 
rate to zero growth at high amplitudes because of the 
limitations of our equipment. 

Our experimental data is well described by a Boltz- 
mann sigmoidal curve for 7 versus F. This fitting rela- 
tion is given by eq. ([!]), and it has four fitting param- 
eters: (a) The higher asymptote jrnax (characteristic 
value for the higher growth rates) (b) The size of the 
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Fig. 3 Growth rate 7 versus the maximum adimensional acceleration F. The dashed hues correspond to sigmoidal fits of each 
curve given by equation (1). Each point represents the average of ten similar experiments, the error bars have been omitted to ease 
visualization. The three arrows in the graph are pointing to the values (T"', 7) that mark an abrupt transition from nonzero 7 to zero 
growth rate at sufficiently large F. 



interval A^y (gives the range of growth rate from the 
lower asymptote to upper asymptote); (c) the inflec- 
tion point Fi] and (d) the width AF of the transition 
region. 

7 = Imax - ^ ^ ^{r-ri)/Ar ' 

A plot of 7max as a function of the frequency is 
shown in Fig. HI This parameter is the saturation growth 
rate and gives an upper limit for the speed of the up- 
ward granular transport through the tube. We observe 
that the ability of the granular material to flow up- 
wards is reduced as the frequency increases. For higher 
frequencies we can expect that the granular system ap- 
proaches a liquid like behavior, and consequently, as the 
frequency rises it becomes more difficult for the granu- 
lar bed to climb up. 

For each frequency we can separate the observed be- 
havior in Fig.[3]in three regions. One for which 7 = at 
low F^ the second one for which 7 7^ 0, for intermediate 
values of F] and a third region for which 7 = and F 
is high. We can speculate that, for a given frequency, at 
low F] the granular bed behaves like a solid; at inter- 
mediate values of F] a mixed state occurs were the sohd 
like state and the fluid state combine in a complex way 
to produce the collective upward motion of the grains 
in one of the branches of the container; and for suffi- 
ciently high F the liquid like state dominates and the 
grains are unable to climb up on one of the branches of 
the tube. 
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Fig. 4 Plot of the maximum growth rate ^max as a function of 
frequency. It is clear that higher frequencies reduce the ability of 
the granular material to flow upwards. 



We will now focus below on two types of transitions 
that were observed. One of them was a (seemingly con- 
tinuous) transition from zero growth at the lowest am- 
plitudes to a finite 7 at higher amplitudes. The other 
was a "jump" in the value of 7 (resembling a first order 
phase transition), which marked a transition from finite 
growth rate to no growth, when F was sufficiently large. 
A detailed study of each transition requires further in- 
vestigations. For the case of the abrupt transition, a 
different experimental setup capable of reaching larger 
amplitudes is necessary. Nevertheless, our data could 
give us an insight on the nature of the transition at low 
F. 
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To probe the nature of the transition at low we 
plot in Fig. O r = -fmaxh versus ^ = (r - rc)/rc. 
The parameter r can be thought as a characteristic re- 
laxation time for the evolution from a symmetric to an 
asymmetric configuration. The critical value Fc of the 
adimensional acceleration was taken as the mean be- 
tween the highest F before 7 becomes non zero, and 
the lowest value of F for which 7 is different from zero. 
At the lower amplitudes it was observed that the critical 
value Fc for the transition from zero growth to nonzero 
7, shifted to higher values of F as the frequency in- 
creased. We can see how r diverges as F approaches 
the critical value Fc. More experiments have to be done 

be^"" .^ „r.u:„ 




Fig. 5 The U-tube instability seems to exhibit a continuous 
phase transition. Here we can see how the characteristic time 
r(^) diverges as F approaches a critical value Fc ~ 0. The data 
shown correspond to all the non zero 7 measured. 



to the stable fixed point). The smaller 1/AF^ the more 
energy needs to be supplied through vibration to reach 
the plateau of the curve 7n(^)- In Fig. [6] we plotted 
the width AF of the region where 7 grows. We observe 
that AF rises with the frequency up to / = 12.5 Hz 
and then drops almost an order of magnitude. This be- 
havior signals a qualitative change around / = 12.5 Hz. 
for this frequency 1/AF is much smaller than for the 
other ones we measured, therefore around this partic- 
ular frequency we need to supply more energy to the 
granular bed to make it go from zero growth rate to 
the fastest growth. 
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Fig. 6 Plot of AF versus the frequency of vibration. The param- 
eter AF is a measure of the increase in energy supplied to the 
granular bed through vibrations, to go from zero growth to the 
plateau of the curve j{F) (see fits in Fig.[3l) for a given frequency. 
A sharp increase occurs around 12.5 Hz. 



Finally, we want to concentrate on the region of 
nonzero growth, for each curve we fixed the frequency 
and changed the amplitude A of vibration, assuming 
that 7min ^ 0, we can rewrite equation ([T]) in the fol- 
lowing form: 



7n(r) = 



1 



1 



o{r,-r)/Ar' 



(2) 



where 7n 1 hmax- 

This expression gives the value of the normalized 
growth rate as a function of F . For each frequency 
equation (j2j) is a solution of the following differential 
equation: 



dF AF ^ 



■7nj7n. 



(3) 



This equation is analogous to the well know logis- 
tic equation so we can picture 7^ = 1 as a stable fixed 
point for jri{F). Therefore as F increases the growth 
rate approaches its maximum value (which correspond 



A physical model that would explain our data and 
the above relationship for the growth rate and the adi- 
mensional acceleration for a given frequency is beyond 
the scope of this paper, but we consider that the sim- 
plicity of the results should help to suggest the rele- 
vant conditions that need to be satisfied by a reasonable 
model for the behavior of a vertically shaken granular 
bed in a partially filled U-tube. 



4 Conclusions 

We have made an experimental characterization of the 
rate of growth for the height difference Ah between 
the free levels of the granular columns in a partially 
filled U-tube subject to vertical vibrations. We have 
observed a continuous transition from zero growth to 
non zero growth rate 7. The region of finite growth 
was well described by a differential equation analogous 
to the logistic equation. We also observed an abrupt 



transition from high values of 7 to zero growth similar 
to a first order phase transition. 
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